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The four-body bound state with two-body forces is formulated by the Three-Dimensional 
approach, which greatly simplifies the numerical calculations of few-body systems with- 
out performing the Partial Wave components. We have obtained the Yakubovsky equa- 
tions directly as three dimensional integral equations. 



1. Introduction 

The four-body bound state calculations are traditionally carried out by solving 
coupled Yakubovsky equations in a partial wave basis. After truncation this leads 
to two coupled sets of finite number of coupled equations in three variables for 
the amplitudes. This is performed in configuration space[l] and in momentum 
space[2,3]. Though a few partial waves often provide qualitative insight, modern 
four-body calculations need 1572 or more different spin, isospin and angular mo- 
mentum combinations [3]. It appears therefore natural to avoid a partial wave rep- 
resentation completely and work directly with vector variables. This is a common 
practice in four-body bound state calculations based on other techniques [4-9]. In 
recent years W. Glockle and collaborators have introduced the three-dimensional 
approach which greatly simplifies the two- and three-body scattering and bound 
state calculations without using partial wave decomposition[10-12]. In this paper 
we extend this approach for four-body bound state with two-body interactions, we 
work directly with momentum vector variables in the Yakubovsky scheme. As a 
simplification we neglect spin and isospin degrees of freedom and treat four-boson 
bound state. Although the four-boson bound state has been studied with short- 
range forces and large scattering length at leading order in an effective quantum 
mechanics approach[13], but it is also based on partial wave approach. 

2. Momentum Space Representation of Yakubovsky Equations in 
3-D approach 

The bound state of four identical bosons which interact via pairwise forces is given 
by coupled Yakubovsky equations [13]: 
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|V'l> = G il2P[(l + P34)|^l) + |^2)] 

\*h) = G tl2P[{l + PMl) + Ifa}] (1) 

In order to solve coupled equations, Eq.(l), in momentum space we introduce the 
four-body basis states corresponding to each Yakubovsky component: 

|ui u 2 u 3 ) 

W1V2V3) (2) 

Let us now represent coupled equations, Eq.(l), with respect to the basis states 
have been introduced in Eq.(2): 

(Hi u 2 «3|V'i} = J D 3 u J D 3 u (ui u 2 uz\GotP\u 1 u 2 u 3 ) 

(u[ u 2 u 3 \(l + P 34 )|ui u 2 u 3 ) {u[ u 2 M3 l^i) 
+ J D 3 u J D 3 v (uiU2U3\GotP\u[u 2 u 3 ) 

{uxu 2 v 2 v 3 ) (t?!^ #3 1 V'2) 
{v\v 2 vz\ip 2 ) = j D 3 v j D 3 u (vi^vslGotPl^^v^ 

(vi v 2 v 3 \(l + P 3i )\u 1 u 2 u 3 )(u 1 €t 2 U3IV1) 

+ J D3v ' («1^2«3|G tP|w 1 l/ 2 4}«l44}(«l4«3lV'2) (3) 

Where D 3 A = d 3 Ai (PA2 d 3 A%. After evaluating the following matrix elements: 



(Ul U2 uz\G$tP\u 1 u 2 u 3 ) 



S 3 (u 3 -u 3 ) 



m 4m 3m 



{ S 3 (u 2 + u 1 + ^u 2 ) (ui|i(e)|^u 2 + u 2 ) 
+S 3 (u 2 - u[ + iu 2 ) (wi|i(e)|-^u 2 - u 2 )} 

JTJ _ Zl 12_ _ ^3 

m 2m m 

(u 1 u 2 u 3 \(l + P 3 4)|mi u 2 u 3 ) = S 3 (u 1 - w'i) 

x { 5 3 (u 2 - u 2 )S 3 (u 3 - u 3 ) 

/ 1 8 / 1 

$ 3 (u2 - 3^2 - gU 3 ) 8 3 (u 3 -u 2 + -u 3 ) } 

{v\V 2 v 3 \{l + P 34 )\u 1 u 2 u 3 ) = S 3 ^ -vj 



1 2 1 2 1 1 1 / / 

x { 8 3 (u 2 + -v 2 - -v 3 ) 5 3 (u 3 + -v 2 + v 3 ) 

+S 3 (v! 2 + -v' 2 + -v' 3 ) 5 3 (u' 3 + -v' 2 - v' 3 ) } (4) 
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We can rewrite the coupled equations, Eq.(3), as below coupled three- 
dimensional integral equations: 

1 f i 1 i 

(ui^wsl^i) = ~ x ^2 / d 3 u 2 {u 1 \t s {e)\-u 2 + u 2 ) 

m 4m 3m 

x { ("2 + t;u 2 u 2 u 3 |^i) 

1 1 _/ 8^ v 1 _ , , . 

+ (U 2 + -M 2 -M 2 + -ti 3 M 2 - jj«3|W 

+ (U2 + 2 M 2 - «2 - 2 M 2 - 3^31^2) } 



{viv 2 v 3 \ip2} = - 2 2 v2 ~r I rf 3 u 3 (tTi|i s (e*)|4) 

m 2m m 

2^ 2/1 



x { 2 (v 3 -v 2 + -v 3 -v 2 - v^i) 

+(V3 - v 2 vz\if) 2 ) } (5) 

Here e and e* are two-body subsystem energies and (a\t s (e)\b) is the sym- 
metrized two-body i-matrix[10]. The so obtained Y-amplitudcs fulfill the below 
symmetry relations, as can be seen from Eq.(5). 



(uiU 2 U3\i>l) = (-UlU2U3\i>l) 
(Hi V 2 V3\lp 2 ) = (-Vi V 2 V 3 \lp2) 

(Vl V 2 V 3 \lp 2 ) = {Vi V 2 - V3\^ 2 ) (6) 



3. Choosing Coordinate Systems 

The Y-components \ipi(a b c)) are given as function of Jacobi momenta vectors as 
solution of coupled three-dimensional integral equations, Eq.(5). Since we ignore 
spin and isospin dependencies, for the ground state both Y-components \ipi(a b c)) 
are scalars and thus only depend on the magnitudes of Jacobi momenta and the 
angles between them. The first important step for an explicit calculation is the 
selection of independent variables, one needs six independent variables to uniquely 
specify the geometry of the three vectors[12]. Therefore in order to solve Eq.(5) 
directly without introducing partial wave projection, we have to define suitable 
coordinate systems. For both Y-components we choose the third vector parallel to 
Z— axis, the second vector in the X — Z plane and express the remaining vectors, 
the first as well as the integration vectors, with respect to them. With this choice of 
variables we can obtain the explicit representation for the Y-components \ipi) and 
|V 2 ) as [14]: 
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ip 2 (vi v 2 v 3 X 1 X 2 Xf 2 ) 



Tpi(u 1 u 2 U 3 XlX 2 x\ 2 ) 



x 




x 



{ 2Mv 3 £i S 2 X 12 X 13 X^J 



ip2{v 3 v 2 v' 3 X 22 X 23 Xv%v 2 ) } 



(7) 



The above coupled three-dimensional integral equations are the starting point for 
numerical calculations. 

4. Summary 

An alternative approach for four-body bound state calculations, which are based 
on solving the coupled Y-equations in a partial wave basis, is to work directly with 
momentum vector variables. We formulate the coupled Y-equations for identical 
particles as function of vector Jacobi momenta, specifically the magnitudes of the 
momenta and the angles between them. We expect that coupled three-dimensional 
Y-equations can be handled in a straightforward and numerically reliable fashion. 
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